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ABSTRACT 

This  report  reviev/s  the  possibility  of  large  airblast 
attenuation  in  the  case  of  rough-walled  tunnel  systems 
in  underground  ammunition  storage  sites.  A comparison  is 
made  betv/een  the  peak  pressure  attenuation  observed  in  three 
tests  with  tunnel  diameters  ranging  from  2,7  - 6m.  These 
results  are  shovm  to  be  consistent  with  a semi-empirical 
theory  for  steady  turbulent  flov;  in  rough  tunnels. 

The  use  of  tunnels  with  large  wall  roughness  is  shown  to 
reduce  considerably  the  hazardous  area  around  an  vinderground 
ammunition  storage  site  compared  to  sites  v/ith  smooth-walled 
tunnels . 
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1 . INTRODUCTION 

1 . 1 Background 

This  report  originated  in  connection  with  the  analysis 
of  a large  scale  test  called  "Operation  Block"  /1,  2/  ■ 
concerning  safety  aspects  of  underground  ammunition  stor- 
age. One  part  of  this  test  was  to  measure  the  blast  at 
various  distances  along  the  main  passagev;ay  leading  into 
the  storage  chamber.  These  results  shov;ed  such  a high 
shock  front  attenuation  compared  to  the  results  from 
model  tests  /3,  4/  that  it  initially  caused  a suspision 
that  there  had  been  erroneous  readings  from  the  pressure 
gauges.  However,  subsequent  analysis  shov^ed  that  the 
results  v^ere  qualitatively  consistent  with  the  large 
difference  in  wall  roughness  between  the  prototype  and 
the  model  /3,  4/. 

1 . 2 Objectives 

The  basic  objective  of  this  report  is  to  show  that  the 
large  airblast  attenuation  observed  in  the  "Operation 
Block"  test,  is  in  fact  comparable  to  the  results  from 
various  other  large  scale  tests  v/ith  similar  wall  rough- 
ness. These  data  are  thus  of  considerable  interest  as 
there  are  strong  indications  that  the  hazardous  area 
around  an  underground  ammunition  storage  site  will  be 
reduced  considerably  if  the  tunnel  walls  have  a large 
roughness. 

1 .3  Scope 

In  Sec.  2,  a phenomenological  description  of  airblast  in 
tunnels  is  given  and  a semi-empirical  theory  for  steady 
turbulent  flov/  in  rough  tunnels  is  reviev/ed.  This  v/ill 
be  compared  with  experimental  findings  in  Sec. 3.  The 
effects  of  tunnel  roughness  on  external  safety  distances 
in  underground  ammunition  storage  sites  will  be  discussed 
in  Sec. 4 and  the  principal  results  are  finally  summarized 
in  Sec. 5. 
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2.  THEORETICAL  CONSIDERATIONS 

A great  deal  of  attention  has  been  directed  in  recent 
years  to  the  subject  of  shock  front  attenuation  due  to 
wall  roughness  as  evidenced  by  the  experimental  and 
theoretical  work  reviewed  in  Ref. 5.  For  the  most  part, 
these  tests  have  been  performed  in  relatively/  smooth- 
walled  tubes  (shock  tubes)  of  moderate  diameters  (order 
of  centimeters).  There  are,  hov/ever,  various  data  on  air- 
blast  in  rock-tunnels  of  relatively  large  diameters  (order 
of  meters)  with  considerable  wall  roughness.  These  data 
are  therefore  of  interest  in  connection  with  underground 
ammunition  storage, 

• Before  specific  examples  are  discussed,  a brief  reviev/ 
will  be  given  of  some  of  the  effects  influencing  the 
blast  wave  propagation  in  straight  tunnels  such  as: 

a)  The  friction  between  the  moving  air  and  the  tunnel 
wall, 

b)  The  pressure  gradients  in  the  wave  behind  the  shock 
front  (rarefaction  v/aves)  v/hich  continuously^’  degrades 
the  front, 

c)  . Thermal  energy  transmission  to  the  walls  (conduction/ 

convection  and  radiation). 

d)  Elastic  or  inelastic  deformation  of  the  walls. 

For  rough-walled  rock  tunnels,  the  blast  wave  attenuation 
is  dominated  by  a) . For  peaked  v/ave  forms  the  effect  in 
b)  will  also  be  of  some  importance.  The  effects  in  c)  and 
d)  are  expected  to  be  comparably  imimportant  as  discussed 
elsewhere  /6/. 

In  the  follov/ing  sections,  various  aspects  of  the  initial 
blast  wave  formation  in  tunnels  will  be  discussed 
follov/ed  by  a discussion  of  the  effects  from  a),  assuming 
a relatively  flat  pressure  wave.  In  this  case  the  effect 
in  b)  can  be  neglected,  and  will  be  discussed  separately 
in  Sec.  2.6. 


k. — . 
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2.1  Blast  V/ave  Formation  in  T^annels 

^ ,^M  

Before  a steady  state  one-dimensional  blast  v/ave  propaga- 
tion in  tunnels  takes  place,  there  are  various  initial 
blast  wave  diffraction  patterns  characteristic  for  the 
specific  geometries  in  question.  Typical  examples  are 
shov/n  schematically  in  Pig.  2.1  with: 

a)  Detonation  of  a charge  in  a straight  tunnel. 

b)  Shock  v;ave  entering  a tunnel. 

c)  Detonation  of  a charge  in  a chamber  adjacent  to  the 
tunnel . 

Cases  a)  and  c)  are  of  particular  interest  for  under- 
ground ammunition  storage,  whereas  case  b)  is  of  interest 
in  underground  protective  installations. 

Porzel  111  has  given  a phenomenological  description  of 
case  b)  in  Fig.  2.1  v/ith  a surface  v/ave  entering  a tunnel 
and  suggested  that  the  subsequent  blast  wave  propagation 
down  a rough-walled  tunnel  is  characterized  by  five 
mechanisms: 

Diffraction,  buildup,  turbulent  choke,  wall  friction  and 
feedback. 

t» 

Although  the  situation  is  somewhat  different  for  cases 
like  a)  and  c)  relevant  for  underground  ammxmition  stoi'age, 
this  model  description  will  give  some  insight  into  the 
basic  phenomena  and  will  therefore  be  reviewed  briefly  in 
the  following  sections. 

2.2  Diffraction  Region 

The  initial  blast  wave  propagation  in  the  timnel  is 
characterized  by  complicated  diffraction  patterns.  It 
is  therefore  not  readily  possible  to  define  a unique 
diffracted  pressure  p^^.  For  engineering  use  it  is  at 

best  possible  to  define  p^j  as  an  average  pressure  in 
space  and  time. 


■ t . 
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2.3  Build-up  Region 

After  a fev/  tunnel  diameters  (l/D)^,  say,  there  is  a 

pressure  huild-up  due  to  the  succession  of  reflected 
pressure  waves  off  the  tunnel  v/alls  v/ith  the  formation  of 
a Mach  stem.  Again,  it  is  hardly  possible  to  define 
a build-up  pressure  p^  at  a specific  distance  (L/D)^. 

For • engineering  use  p^  can  approximately  be  thought  of  as 
being  the  average  pressure  of  more  or  less  definite  shocks 
v/hich  coalesce  at  a fev/  tunnel  diameters  dov/nstreara  from 
the  diffraction  region.  For  example,  for  case  b)  in 
Fig.  2.1  with  a surface  v/ave  entering  a tunnel,  empirical 
fits  to  various  measurements  have  tentatively  produced 
Pg  Pi  (1,4  i 0,2)  pjj  /?/. 

To  calculate  these  types  of  diffraction  patterns  from 
first  principles  for  a real  three-dimensional  case,  appear 
to  be  prohibitively  difficult.  It  is  therefore  usually 
necessary  to  rely  on  previous  experiments  or  specific 
model  tests  in  each  particular  case. 


2.4 


2.5 


Turbulent  Choke 

After  the  formation  of  the  more  or  less  definite  Mach  stem 
described  in  Sec.  2.3,  the  wall  roughness  produces  turbu- 
lent boundary  layers  which  steadily  -grow  out  from  the  v/all. 
Soon  the  boxindary  layers  will  converge  at  a certain  dis- 
tance behind  the  shock  as  indicated  in  Fig.  2.4.  At  the 
center  of  this  closure  the  material  flow  v/ill  be  stagnated. 
This  situation  is  refered  to  as  a turbulent  choke.  Assum- 
ing a steadj'-- state  situation  with  a constant  distance  for 
choke  formation  behind  the  shock  and  isentropic  flow  be- 
hind the  shock,  the  model  relates  the  "choked"  overpressure 
p^  at  the  front  to  the  buildup  overpressure  Pg  as 


Pb  = (1  + Po) 


1 + 


7/2 


- 1 


(2.4) 


with  pressures  expressed  in  bars. 

Shock-Front  Attenuation  due  to  b-all  Roughness 
An  example  of  a shock  v/ave  passing  over  an  idealized 
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roughness  element  is  shorn  in  Pig.  2.5a.  Clearly,  in  a 
naturally  excavated  tunnel  v;ith  large  wall  projections 
at  random,  there  v/ill  he  a very  complex  pattern  of  shock 
reflections  which  v;ill  drastically  affect  the  shock  front. 
Detailed  calculations  of  the  blast  w'ave  propagation  from 
first  principles  will  obviously  be  prohibitively  difficult. 

It  is  beyond  the  scope  of  this  report  to  give  a critical 
reviev/  of  the  many  approximate  models  dealing  with  this 
problem  /5/.  A so  called  shock  impedance  model  originallj^ 
proposed  by  Porzel  111  is  conceptually  simple  and  appear 
to  reproduce  various  experimental  results  of  the  present 
type  reasonably  well  and  will  be  discussed  next. 

The  basic  assumption  for  this  model  can  be  expressed  as 
f ollov/s : 

The  loss  in  the  shock  energy  due  to  wall  roughness  is 
proportional  to; 

a)  The  kinetic  energy  per  unit  volume  at  the  shock  front, 

TT* 

\,v* 

b)  The  volume  subtended  by  the  average  wall  roughness 
of  the  tunnel  walls  or  ^ • S • dL,  where  e"  = average 
height  of  the  projections,  S = perimeter  of  the 
tunnel,  and  dL  = distance  traversed  by  the  shock. 


c)  An  efficiency  factor  £ postulated  to  be  in  the  range 

1/2  ^ 1. 

The  infinitesimal  loss  in  total  energy  over  a distance 
dL  can  therefore  be  expressed  as 

dE  = - v'^  * ® (2.5a) 

Letting  dS.  _ represent  the  total  energy  per  unit 

X,  V 

volume,  i.e.  dE=  A*  D*dS.  , Ea  (2.5a)  can  be  ex- 

X,  V 

pressed  as 


(2.5b) 


= -4€  (e/D)d(I/D) 

Here,  D designates  the  hydraulic  diameter  defined  by 

D = 

v/ith  A = cross-sectional  area  and  S the  perimeter 
of  the  tunnel . 

Eq  (2.5b)  as  it  stands  is  of  course  not  very  con- 
vinient  for  practical  use,  but  can  be  re-expressed  in 
terms  of  the  front  pressure  using  the  usual  Rankine- 
Hugoniot  equations.  This  produces  an  equation  of  form: 

Y(p)  = constant  - 2 £ (e/D) (L/D)  (2.5d) 

where  Y(p)  is  an  impedance  function  defined  as 

Y(p)  = Inp  - -2^  ln(p  1)  - 1 (2.5e) 

Here,  p is  the  dimensionless  overpressure  ratio 
P = Psq/^o  Pso  ” static  overpressure  and 

p^  = ambient  pressure. 

Expressing  p^^  in  bars  and  using  p^  = 1 bar,  we  have 
nominally  p = P^q*  V/e  will  therefore  in  the  remainder 
of  this  report  not  distinguish  between  p and  p_*. 

Furthermore,^^  is  a thermodynamic  parameter 

i)/(  y - ^)  (2.5f) 

where  y is  the  usual  ratio  of  the  specific  heats, 
y = Cp/C^.  For  air  under  standard  conditions  y = 1,4- 

which  gives  yCo=  6 and  Eq  (2.5e)  takes  the  simpler 
form 

Y(p)  = ^^-Inp  - ^ln(p  + 7)  - ^ . (2.5g) 

This  function  is  shown  in  Fig  2.5b. 


L. 
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It  is  interesting  to  note  that  in  the  high 
pressure  limit,  p»  7,  Eq  (2.5g)  reduces  to 

Y(p)  = Inp  and  therefore  Inp  = constant  - 2£  ^ , 


which  has  the  solution: 


?2  = P.|exp  (-k  g ), 


(2.5h) 


v/ith  an  attenuation  constant 


(2.5i) 


Eq  (2.5h)  is  the  form  for  the  blast  wave  attenuation 
originally  proposed  by  Emerich  and  Curtis  in  Ref, 3, 
This  equation  was  also  used  by  this  office  in  an 
earlier  analysis  of  blast  v;ave  propagation  in  tubes 
and  tunnels  /9/. 

Eq  (2.5d)  is  in  a form  which  allows  the  determination 
of  pressure  attenuation  due  to  a given  relative 
roughness  e/S,  The  length  of  tunnel  required  to 
reduce  the  pressure  from  p^  to  P2  can  thus  be  found 
from 


This  shows  that  there  ought  to  be  a linear  rela- 
tionship betv/een  the  impedance  function  Y and  the 
distance  L in  the  tunnel. 

The  definitions  of  e and  £ need  futher  clarifica- 
tions. For  very  irregular  tunnel  walls,  the  deten 
mination  of  e will  obviously  be  quite  uncertain. 
This  is  of  course  not  only  the  case  for  this 


Y(p^)  = constant  ~2£^  ^ 


(2.5j) 


Y(P2)  = constant  “26  § ^ 


(2.5k) 


(2.51) 
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particular  model,  "but  for  any  model  attempting  to  describe 
wall  roughness  attenuation.  Por  engineering  use,  an 
adequate  value  can  be  loosly  expressed  as 


where  e equals  the  depth  of  the  roughness  elements 
and  the  area  is  integrated  over  the  wall  of  the 
tunnel . 

Examples  of  idealized  wall  roughness  elements  are 
shovai  in  Pig.  2. 5c  and  estimates  of  e using  Eq  (2.5m) 
are  given.  Por  a tunnel  blasted  out  of  rock,  a 
reasonable  value  for  e"  appears  to  be 


vdiere  e in  the  maximum  height  of  the  roughness 
elements.  Clearly,  only  experiments  v/ill  be  able  to 
clarify  the  adequateness  of  this  definition  of  e. 

This  v/ill  be  discussed  in  Sec. 3. 

The  constant  ^ in  Eq  (2.5d)  is  related  to  the 
amount  of  kinetic  energy  stagnated  within  the 
volume  extended  by  the  roughness.  It  has  been 
postulated  that  £ = 1 in  the  early  period  of  the 
shock  propagation  of  the  tunnel,  i.e.  all  the 
kinetic  energy  is  stagnated  and  not  immediately 
available  for  supporting  t*he  shock  ("strong  impedance"). 
At  later  times,  approximately  50/b  of  this  stagnated 
energy  becomes  hydrodynamically  avaliable  to  be  pro- 
pagated in  the  forv/ard  direction  toward  the  shock 
front. 

This  is  therefore  a "feedback"  situation  and  one 
has  approximately  € = (’’weak  impedance").  Earlier 
tests  have  indicated  that  there  is  a more  or  less 
definite  distance  at  which  point  £ changes  its  value 
from  1 to  . Empirically,  this  "feedback"  distance 


(2.5a) 
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has  been  determined  to  be  approximately  given  by: 

_ 

(L/D)-,  18/(6  |) 


(2.5n) 


For  an  effective  roughness  factor  0,05  e/D<0,1, 

the  corresponding  feedback  distance  is  44  > (L/D).,i>36. 

2.6  Shock-Front  Attenuation  due  to  Rarefaction  V/aves 

In  the  discussion  of  the  blast  v;ave  propagation  in 
the  previous  sections,  the  input  was  assumed  to  be 
relatively  constant  in  time.  A characteristic 
feature  of  the  propagation  of  a peaked  blast  v;a,ve 
along  a tunnel  is  that  the  higher  pressure  regions 
toward  its  front  move  at  higher  velocities  than  the 
shock  front  and  lov/er  pressure  regions.  The  shock 
front  is  thus  overtaken  by  the  lower  pressures  of 
the  rarefaction  v;ave,  v/hich  causes  an  attenuation 
of  the  shock  front.  At  the  same  time  this  also 
results  in  an  increase  in  the  positive  duration  or 
wavelength  with  travel  distance.  This  is  illustrated 
in  Fig.  2.6a. 

Various  simplified  models  have  been  proposed  to 
characterize  this  phenomenon  /5/,  but  one  of  the 
most  rigorous  descriptions  has  probably  been  pre- 
sented by  Porzel  /?/  with  the  following  result  for 
the  infinitesmal  shock  front  attenuation  with  dis- 
tance: * 


wave  at  the  front  as  indicated  in  Fig  2.6a.  Clearly, 
for  a flat  shock  front,  as  is  approximately  the  case 
in  an  ordinary  air-driven  shock  tube,  ^p/3  t = 0 and 
there  is  no  front  pressure  reduction  due  to  this 


(2.6a) 


Here,  f(p)  is  a relatively  complicated  function 
of  the  front  pressure  p and  u is  the  shock  front 
velocity,  ^p/  c)t  is  the  slope  of  the  rarefaction 


effect. 


The  use  of  Eq  (2.6a)  for  peaked  shock  waves  presumes 
a knowledge  of  the  initial  pressure-time  history  at 
some  position  in  the  tunnel.  In  practice,  for  under- 
ground ammunition  storage,  Bp/2)t  cannot  easily  he 
defined.  A purely  empirical  approach  based  on  model 
tests,  has  therefore  been  used  to  establish  the  "smooth  . 
v/all"  pressure-distance  relationships  for  various  explo- 
sive weights  in  a relatively  wide  range  of  tunnel 
systems  shown  in 'Fig  2.6b  /10-12/: 

^2 

p=  f(a^,  a2»  ...)(Q/Vj.)  ^ (2.6b) 

Here,  Q is  the  explosive  charge  (TNT)  and  the 
total  volume  traversed  by  the  blast  wave  up  to  the 
observation  point  as  indicated  in  Fig  2.6b. 
f(a^,  ag,  ...)  in  Eq  (2.6b)  is  a function  depending 
on  the  exact  geometry  around  the  area  of  detonation 
with  a^ , ag,  ...  indicating  geometrical  parameters 
such  as  cross-sectional  areas,  angles  etc.  The 
exponent  0^  is  found  to  be  remarkably  reproducible 
for  the  various  configurations  with  O2  = “ Of"!. 

Eq  .(2. 6b)  shows  therefore  that  for  a given  configura- 
tion the  pressure  varies  roughly  as 

po=  (Q/V^)°»^  ■ . (2.6c) 

These  observations  will  be  used  later  in  Sec  3 
in  the  analysis  of  blast  wave  attenuation  in  rough- 
walled  tunnels. 

It  is  interesting  to  compare  the  pressure  dependence 
for  various  confinements.  In  Fig  2.6c  the  results 
for  TNT  detonations  in  free  air  /14/  in  a tunnel  /IO/, 
and  in  a closed  chamber  /13/  are  shown.  Here,  the 
average  peak  pressures  are  given  versus  an  effective 
loading  density  Q/V,  where  Q is  the  TNT  charge  v;eight 
and  V the  total  volume  traversed  by  the  blast  wave 
up  to  the  observation  point.  As  may  be  seen,  the 
results  differ  remarkably  little  in  the  pressure 
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range  5-100  bar  corresponding  to  effective  loading 
densities  ranging  from  0,1  to  10  kg/m3.  In  fact, 
the  differences  are  less  than  about  50vb,  For  higher 
and  lower  loading  densities  the  differences  are  be- 
coming significant.  It  is  also  interesting  to  note 
that  the  tunnel  data  on  the  average  appear  to  lay 
roughly  between  the  free  air  and  closed  chamber  data. 
This  is  reasonable  on  physical  grounds  because  TNT 
is  deficient  in  oxygen  and  the  detonation  energy  is 
dependant  on  the  confinement. 


3.  EXPERIMENTAL  RESULTS  AI^D  DISCUSSION 

3 . 1 Large  Scale  Tests  in  Tunnels 

A full  scale  test  called  ’’Operation  Block”  v/as  con- 
ducted in  1973  by  the  Royal  Swedish  Fortification 
Administration  /l,  2/.  The  basic  objective  of 
this  test  'was  to  study  the  possibility  of  reducing 
the  hazardous  area  outside  an  underground  ammuni- 
. tion  storage  by  installing  a moveable  concrete 
block  in  the  branch  passageway.  In  conjunction  with 
the  full  scale  test,  a series  of  model  tests  in  a 
linear  scale  1:50  were  performed  in  Norway  /3,  4/. 

In  particular,  one  part  of  the  tests  was  to  measure 

the  pressure-time  history  in  the  branch-  and  main  passageway 

leading  into  the  storage  chamber  as  shown  in  Fig  3.1a. 

The  results  of  these  measu’rements  are  shovm  in  Fig  3.1b. 

As  may  be  seen,  the  model  tests  showed  a relatively 
small  shock  front  attenuation  in  the  straight  tunnel 
section  (MP  7-9)  whereas  the  full  scale  results 
showed  a considerable  reduction.  This  imexpected 
result  caused  initially  a suspicion  that  the  reading 
from  the  pressure  gauge  nearest  to  the  exit  had  been 
erroneo’as. 

Hov/ever,  subsequent  calculations  of  the  front 
pressures  from  arrival  time  data  confirmed  to 
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within  approximately  i 20/o  the  earlier  average 
peak  pressure  readings  /15/  as  seen  from  Fig  5.1b. 

The  large  attenuation  observed  in  the  outer  parts 
of  the  tunnel  has  been  analysed  elsewhere  /3/  and 
found  to  be  qualitatively  due  to  the  large  v/all 
roughness.  This  will  be  discussed  shortly,  but 
first  we  v/ill  show  the  results  from  two  other  large 
scale  tests  in  tunnel  systems  with  similar  wall 
roughness.  The  geometries  and  the  positions  of  the 
measuring  points  in  these  tests  are  shov/n  in 
Pig  3.1c  and  3 . 1 d . 

The  tunnels  were  blasted  out  of  rock  with  hydraulic 
diameters  approximately  D = 2,7  m /16/  and  D = 6 m 
/IO/,  The  average  height  of  wall  roughness  elements 
v/as  approximately  e = 0,2  - 0,1  m (defined  in 
Sec  2.5)  and  comparable  to  the  v/all  roughness 
in  the  "Operation  Block"  test. 

In  the  2,7  m diameter  tunnel  system.,  TNT  charges 
ranging  from  100  kg  to  1000  kg  v/ere  detonated  near 
the  dead  end  of  the  tunnel,  v/hereas  in  the  6 m dia- 
meter tunnel  TNT  charges  ranging  from  1 kg  to  92  kg 
were  detonated  near  one  of  the  exits.  Fig  5.1e 
shows  typical  peak  pressure  data  versus  distance  in 
tunnel  diameters  for  these  tests (defined  in  Figs  3.1a, 
3.1c,  and  3.  Id)  together  v/ith  the  "Operation  Block" 
data. 

For  details  regarding  data  reduction,  reference  is 
made  to  the  original  reports.  It  should,  hov/ever, 
be  pointed  out  that  the  average  peak  pressure  data 
in  Fig  3.1e  were  found  directly  from  the  pressure- 
time recordings  by  taking  the  mean  values  of  the 
pressures  of  peaks  and  troughs  of  the  pressure 
oscillating  system.  The  data  found  from  shock  front 
velocity  measurements  will  also  correspond  qualita- 
tively to  this  definition. 
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As  may  be  seen  from  Pig  5.1e,  the  blast  attenuation 
is  significant  for  all  three  tests  and  as  v/ill  be 
discussed  shortly,  a large  part  of  this  is  due  to 
a relatively  large  v/all  roughness.  Unfortunately, 
the  results  extend  out  to  moderate  distances  in 
tunnel  diameters  and  this  complicates  a direct 
comparison  with  the  model  discussed  in  Sec  2.  To 
get  further  insight  into  the  long  range  effects  of 
v/all  roughness  attenuation  \ie  shall  first  turn  to 
some  earlier  test  results  in  straight  tubes,  and 
return  in  Sec  5.3  to  the  data  in  Pig  3.1e. 

3 . 2 Tests  in  a 0,03  m Diameter  Tube 

The  long  range  effects  of  v;all  roughness  attenua- 
tion have  been  reported  earlier  for  detonations 
inside  a D = 0,05  m steel  tube  where  measurements 
v;ere  extended  out  to  approximately  L/D  = 700  /10/. 

To  analyse  these  earlier  results,  the  peak  pressure 
data  were  fitted  to  the  theoretical  expression  in 
Pig  (2.5d)  and  the  empirical  results  in  Eq  (2.6b) 
as  follows: 

The  smooth-wall  pressure  was  assumed  to  vary  as 

Cp 

Po  = C^(Q/V^)  ^ • (3.2a) 

and  the  wall  roughness  attenuation  was  assumed  to 
be  governed  by  the  impedance  function 

Y(Pq)  - Y(p)  = 26  U . (3.2b) 

Here,  p^  is  the  fictions  smooth  wall  pressure 
at  a distance  L from  the  charge  and  p the  actual 
measured  pressure.  C^,  Op,  and  2£e'/D  v/ere  treated 
as  fitting  parameters.  To  simplify  the ' calculations 
it  was  assumed  that  €=  l/2  (See  Sec  2.5).  The 
results  of  the  fitting  procedure  are  shown  in 
Pigs  3.2a  and  3.2b.  As  may  be  seen,  measurements 
of  the  attenuation  appear  to  be  consistent  with  the 
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proposed  model.  In  particular,  the  average  v'all 
rougnness  v/as  found  to  he  e = 0,46  *•  0,05  mm  v/hich 
compares  favourably  v/ith  the  direct  estimated  value 
of  e = 0,3  mm  /10/,  The  data  points  in  Pig  3.2a 
appear  to  lay  on  a smooth  falling  curve  with  no 
kinks  as  postulated  for  a constant  pressure  wave 

(Sec  2.5).  In  particular,  the  feedback  distance 
using  Eq  (2.5n)  is  found  to  be  (L/D)p  = 75.  The  re- 
sults in  Fig  3.2a  shov/  that  there  is  no  distinct 
kink  at  approximately  this  distance.  If  a constant 
pressure  v;ave  had  been  assumed,  i.e.  p^  = constant 
in  Eq  (3.2a),  the  use  of  Eq  (2.5d)  vrould  have 
approximately  been  consistent  with  a transition 
from  £,=  1 to  1/2  at  L/D  = 178  and  an  e = 0,6  mm. 

, The  results  of  this  analysis  are  shown  in  Fig  3.2c. 
This  shov/s  the  one  has  to  be  careful  in  attempting 
to  analyse  peaked  pressure  waves  only  in  terms  of 
Eq  (2.5d). ‘ 

3 . 3 Analysis  of  Test  Results  in  Tunnels 

Nov;  v;e  shall  return  to  the  large  scale  results  in 
tunnels  in  Fig  3.1e.  The  data  for  the  D = 2,7  m 
and  D = 6 m tunnels  v;ere  analysed  using  Eq  (3.2a) 
and (3. 2b)  as  for  the  D = 0,05  m tube  and  the  re- 
sults are  shown  in  Fig  3. 2d  in  terms  of  the  impe- 
dance function  and  in  Fig  3.2e  using  the  pressures 
explicitly. 

The  average  wall  roughness  was  found  to  be  e = 0,12  m 
for  the  D = 2,7  m tunnel  and  e = 0,15  m for  the 
D = 6 m tunnel. 

This  v/ould  correspond  to  an  average  height  of  the 
maximum  roughness  elements  e = 2e  -=^0,3  m in 
Eig  2.5c.  This  appears  to  be  quite  consistent  with 
direct  observation  of  the  tunnel  walls. 

As  for  the  "Operation  Block"  results  in  Fig  3. 2d  the 
data  are  really  too  few  in  the  straight  tunnel  to 
do  the  same  kind  of  fitting  procedure.  V/e  have 
therefore  tentatively  shown  the  attenuation  one 
would  expect  for  an  average  roughness  of  e = 0,17  m 
or  2 e e’/D  = 0,1. 
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The  "smooth  wall"  curve  a)v/as  estimated  from 
case  d)  in  Fig  2.61.  As  may  be  seen,  the  attenua- 
tion appears  to  be  somewhat  larger  than  the  pre- 
diction 'of  curve  b) . This  may  possibly  be  due  to 
large  amounts  of  dust  and  v/ater  in  the  tunnel  system 
/3/  or  special  effects  due  to  the  particular  geometry 
in  this  test  /17/. 

The  phenomenological  description  of  blast  in  tunnels  also 
included  a diffraction  region,  a build-up  region,  and 
a turbulent  choke.  This  model  applied  strictly  only 
to  a tunnel  adjacent  to  an  infinite  half  space. 

Clearly,  for  detonation  inside  a tunnel  or  in  a 
chamber  adjacent  to  the  tunnel  as  for  the  tests  re- 
ported here,  the  situation  is  somewhat  different. 

In  particular,  the  turbulent  choke  formation  with  a 
build-up  pressure  p^  and  a "choked"  overpressure  Pq, 
v^ould  according  to  Eq  (2.4)  require  a large  pressure 
drop  out  to  20  tunnel  diameters  or  so.  I'll , The 
present  data  do  not  appear  to  support  this  effect. 

The  observed  reduction  in  pressure  seems  to  be  due 
only  to  diffraction,  impedance  (roughness),  and  rare- 
faction, but  it  is  hardly  possible  to  be  more  specific 
at  this  point. 

The  possible  implications  of  v/all  roughness  attenua- 
tion on  external  safety  distances  has  been  discussed 
briefly  in  Ref  3 in  conneotion  with  the  "Operation 
Block"  test,  but  it  is  of  interest  to  give  a more 
general  discussion.  This  will  be  the  subject  of 
Sec  4. 
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4.  SPFECTS  OP  TUNNEL  ROUGHIIESS  ON  EXTERNAL  SAPSTY  DIST.AITCE3 
Measurements  of  the  blast  propagation  outside  models 
of  underground  storage  sites  /18/  and  shock  tubes  /19, 
20/  have  shov/n  that  the  isobar  distances  from  the 
exit  are  primarily  dependent  on  the  exit  pressure. 

In  particular  it  has  been  proposed  rhat  the  external 
safety  distances  or  50  nbar  isobar  distances  can  be 
approximated  by  /2l/: 


D^/D  = 18 

p 

Sector 

0°  - 

30' 

D4/D  = 15 

p 

Sector 

1 

0 

0 

50' 

D^/D  = 12 

p 0-67 

Sector 

60°  - 

90' 

D^/D  = 8 

p 0-67 

Sector 

90°  - 

120' 

D-,/D  =4,5 

p 0-67 

Sector 

1 

0 

0 

C\J 

I80' 

Here,  D is  the  effective  main  passageway  diameter 
in  meters  and  p^  the  exit  pressure  in  bars.  The 
angles  are  measured  from  the  extended  centerline 
of  the  main  passageway  as  sho\vn  in  Pig  4. a,  and 
D^-D^  are  the  safety  distances  in  meters  from  the 
exit. 

It  is  of  interest  to  see  the  effects  of  wall  rough- 
ness attenuation  on  the  safety  dist.an.ces  using  the 
model  outlined  in  Sec  3.2.  Consider  the  connected 
chamoer  storage  site  in  Pig  4.b  which  is  scaled  up 
from  an  earlier  model  in  a linear  scale  1:100  /18/. 

The  model  tests  showed  the  pressure  in  the  main 
passageway  could  be  fitted  empirically  to  /12/: 

Po  = 4,3  (4.b) 

with  Q and  defined  in  Pig  4.b.  This  result  then 
is  valid  for  an  essentially  smooth-walled  tunnel. 
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In  the  event  of  significant  v/all  roughness,  e/D, 
we  shall  use  the  model  proposed  in  Sec  5.2: 

Y(.p)  = Y(Pq)  - (4.c) 

Here,  p^  is  fictitious  pressure  at  a distance 
L/D  defined  in  Hq  (4.h)  and  p the  pressure  at  the 
same  distance  including  the  v/all  roughness.  The 
efficiency  factor  C is  defined  as 

6=1  for  L/D  ^ (L/D)j. 

6=  1/2  for  L/D  > (L/D)p, 

where  (L/D) -5,  = is  the  feedback  distance  determined 
from  Eq  (2,5n).  The  results  of  these  calcula- 
tions are  shown  in  Eig  4.c  for  values  of  e/D 
ranging  from  0,01  to  0, 1 . 

As  an  example,  we  choose  L/D  = 20  and  use  the 
corresponding  pressures  p as  the  exit  pressures 
Pg  in  Eq  (4. a)  to  estimate  the  corresponding  safety 
distances  v/hich  are  given  in  Table  4. a.  This  is 
displayed  graphically  in  Eig  4.d.  As  may  be  seen, 
even  a moderate  wall  roughness  reduces  the  distances 
significantly. 

It  should  be  cautioned  that  these  results  are  based 
on  the  assumption  of  a steady  state  blast  pi*opaga- 
tion  in  the  t\innel.  As  pointed  out  in  Sec  2,  there 
will  obviously  be  a very  complicated  diffraction- 
region  in  the  first  few  tunnel  diameters  near  the 
exit  of  the  branch  passagev/ay.  Hov/ever,  from  the 
analysis  in  Sec  5 the  present  model  appears  to 
qualitatively  reproduce  the  earlier  results  v/hich 
partly  justifies  the  extension  to  the  more  general 
situation  discussed  here. 
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5.  CONCLUSIONS 

The  main  purpose  of  the  present  repo3?t  has  heen  to  pro- 
vide experimental  results  which  confirm  that  a macro- 
scopic surface  rougliness  in  tunnels  attenuate  front 
pressures  of  blast  waves  significantly.  Thus  v;all  roiigh- 
ness  provides  a degree  of  control  for  reducing  the 
hazardous  area  around  an  underground  ammunition  storage 
site. 

Using  a semi-empirical  theory  for  steady  turbulent  flov; 
in  rough  pipes  in  combination  with  empirical  results 
from  model  tests,  it  has  been  possible  to  reproduce  the 
results  in  tunnels  blasted  out  of  rock. 

Various  empirical  cross  checks  of  the  the  results  in 
tunnels  v/ith  diameters  ranging  from  2,7  m to  6 m have 
produced  an  average  v;all  roughness  of  e = 0,15  - 0,05  m. 

For  a typical  underground  ammunition  storage  this  rough- 
ness could  easily  reduce  the  hazardous  area  by  a factor 
of  two  or  more  compared  to  an  essentially  smooth-v;alled 
site.  These  aspects  should  therefore  be  carefully  con- 
sidered by  the  design  engineer  to  reduce  the  safety 
distances  to  a minimum  at  minimum  cost. 
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Table  4. a Safety  distances  outside  an  underground 

ammumition  storage  site  for  various  timnel 
wall  roughness  as  discussed  in  the  text. 


Fig.2,1  Shock  front  diffraction  patterns  shown  scher.iatically  for: 

a)  fetonation  of  a charge  in  a straight  tunnel. 

b)  Shock  wave  entering  a tionnel. 

c)  Detonation  of  a charge  in  an  undergroxmd 
ammunition  storage  chamber, 

ay 


/ /^ / / / / / / yy  J ^ J 1_J / y / / 

Sioclc  j'ront 

■>i’o 


-/  / / -'/’/TV  / T'/V 


^is»2,4  Turbulent  choke  formation  in  a tunnel, 

Pq  is  the  "chocked"  overpressure  at  the  shockfrort, 
and  Pg  is  the  buildup  pressure. 


rig,2*5a  Shock  wave  passing  an  idealized  roughness  element. 
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Fig,2,?b,  Impedance  function  Y(p) 
versus  overpressure  p. 


as  defined  in  Eq,  (2,5g) 
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Pi^,  2-,  5o  Idealized  examples  of  wall  rougliness  nji*  a)  ,1)) , c) , 

®max  maximum  height  of  roughness  elements 

and  ~ the  average  wall  roughness  estimated  from 
Eq.  5m )as earning  e^<D,  For  a realistic  wall  roughness 
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Fig,2,6aln  a)  the  shook  front  is  attonuated  duo  to  the 
rarefaction  waves  wheras  in  o)  (typical  for  an 
aii'-driven  shock  tube)  no  attenuation  v;ill  take 
place. 
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Fig2.6bThe  results  from  various  model  tests  show  that  the  front 
pressure  in  the  tunnel  systems  of  underground  storage 
sites  with  certain  limitations  can  be  expressed 

empirically  as  p = f(a^,a2, )(Q/V^)^2 

Cg  = 0,6  - 0,1  and  a proportionality  factor  f(a^,a2,,..) 
containing  geometrical  parameters*  ' 
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^g.2,6c  Static  overpressure  versus  effective  loading  density 
for  TNT  detonation  in: 

a)  Free  air  /ih/ • 

b ) Tube  s/ tunnel s /lO/ • 

c)  Closed  chamber  /l3/« 
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Fig.B.lb  j^-verage  peak  pressure  in  a tunnel  system  versus 

distance  from  the  entrance  to  the  detonation  chamber 
in  the  "Operation  Block"  . tests,  seeJ'ig*3*l®. - '-C  fui.l 
scale  data  show  a much  larger  blast  v/ave  attenuation 
near  the  tunnel  exit  than  the  model  data. 


CROSS- SECTONAL  AREA  A = 5,5in^ 
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Tuim©!  system  and  measuring  points  in 
the  Raufoss  trials  fib  /, 
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Blast  propagation  test  in  an  abandoned  railway 
tunnel  ( radius  of  curvature  approximately  200in  ) 
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Fig.3.2a  Peak  pressure  attenuation  of  the  blast  from  9,5g  TNT' 
detonated  in  a 5 cm  diameter  tube  expressed  as  the 
impedance  function  Y versus  tube  length  in  diameters. 
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Fig,3.2b  Peak  pressure  attenuation  of  the  blast  from 
17, 5g  TNT  detonated  in  a 5 cm  diameter  tube 
expressed  as  the  impedance  function  Y versus  tube 
length  in  diameters. 


FiS»3*2cTest  of  proposed  impedance  theory  for  TNT  detonation 
inside  a D = 5 cm  diameter  steel  tube  disrep;arding 
attenuation  due  to  the  rarefaction  wave.  The  straight 
lines  represent  least  square  fits  to  the  data  for 
distances  below  and  above  a feedbaok  distance  (L/D)^^,»!'178 
with  an  efficiency  factor  6=1  and  €=i,  respectively. 

The  four  fitted  values  of  e/D  produce  e/D  = 0,0115^0,0008 
for  all  the  data  , or  an  average  wall  roughness 
S = 0,58  - 0,04fflin;  which  is  about  25  °h  higher  than  the 
value  for  e found  in  Pigs*  ^.2a  and  ^■,2b  where  also  the 
effects  from  the  rarefaction  waveiwere  considered. 


Fig.3.2e  The  same  data  as  in  Fig. 3, 2d  except  for  using 
pressures  explicitly. 


Pig,*+*a  The  definition  of  external  safety  distances. 
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Pig.4,b 


T 

1760  la'^,  Gtorago  chamber  volume 
2 

10  ra  , branch  paeaagev/ay  oroae-sectional  area 
2 

20  ra  , main  pas3age\ifay  cross-cectional  area 

5 m,  effective  diameter  of  main  passageway 

10  m,  length  of  branch  passageway 

88000  kg  Till,  Net  Explosive  Quantity 

50  kg/m^,  loading  density 

effective  loading  density 

V ♦ -^jj*  + 2 LAj^  «»  total  effective  volume 

pressure  at  distance  L/D  from  exit  of  branch 
passageway 

pressure  at  the  exit 

Definitions  of  terns  used  in  the  connected  chamber 
storage  site  discussed  in  the  text 


Fig,4,cPeak  pressure  versus  distance  in  tunnel  diameters  for 

the  connected  chamber  storage  site  discussed  in  the  text 
with  relative  average  wall  roughness  e/D  as  parameter. 
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Pig,4,d  Safety  distances  for  an  imderground  ammunition 
storage  site  with  an  essentially  smooth-walled 
tunnel  in  comparison  with  tunnels  having  a relative 
wall  roughness  e/h  ranging  from  0,02  to  0,1. 
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